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We propose and demonstrate an effective mode- filtering technique of non- Gaussian states gen- 
erated by photon-subtraction. More robust non- Gaussian states have been obtained by removing 
noisy low frequencies from the original mode spectrum. We show that non- Gaussian states preserve 
their non-classicality after quantum teleportation to a higher degree, when they have been mode- 
filtered. This is indicated by a stronger negativity —0.033 =b 0.005 of the Wigner function at the 
origin, compared to —0.018 =b 0.007 for states that have not been mode- filtered. This technique can 
be straightforwardly applied to various kinds of photon-subtraction protocols, and can be a key in- 
gredient in a variety of applications of non- Gaussian states, especially teleportation-based protocols 
towards universal quantum information processing. 



I. INTRODUCTION 

In the field of continuous- variable (CV) quantum in- 
formation processing (QIP), non-Gaussian states of light 
have attracted great interest recently as an essential re- 
quirement for universal QIP, which is unattainable only 
with Gaussian states and operations [1]. To generate a 
variety of non-classical non- Gaussian states, the photon- 
subtraction technique has been playing a central 
role. With this technique, more advanced and powerful 
CV protocols with non-Gaussian states have been pro- 
posed and also experimentally demonstrated recently 

In practice, such CV protocols are less efficient in the 
lower frequencies around the laser carrier frequency due 
to various sources of noises, as well as higher frequen- 
cies due to the finite bandwidths associated with squeez- 
ing, and also the electronics used. To achieve the robust 
implementation of these protocols, fragile non-Gaussian 
states should be encoded properly within the highly- 
efficient frequencies. However, such encoding has been 
a challenging task in photon-subtraction, where the fre- 
quency mode of non-Gaussian states is determined by 
a mode of squeezed light associated with probabilistic 
photon-detection. Sideband- mode encoding (Fig. |l(a)[ ), 
often used for Gaussian states, can also be achieved 
here with non-degenerate squeezed modes, but it re- 
quires high-speed electronics to observe the generated 
states [il. In contras t, cen ter-band encoding with a de- 
generate mode (Fig. |l(b)[ ) is subject to low-frequency 
noises, which need to be removed with electrical high- 
pass-filters (HPF) during homodyne-detection. These 
HPFs also filter out the low-frequency information of 
non- Gaussian states, thereby causing an inevitable fre- 
quency mode mismatch which degrades the fragile non- 
Gaussian states. 



Here we propose a mode-filtering strategy to achieve 
a more efficient encoding of non-Gaussian states via 
photon-subtraction. Our strategy is to filter out op- 
tically the problematic low-frequencies of the degener- 
ate squeezed mode, instead of using its full spectrum 
(Fig. l(c)[ ). By introducing an optical HPF with any de- 
sired cutoff in the photon-detection channel, we can tailor 
the frequency mode of the non- Gaussian states without 
any direct loss on or manipulation of the state itself. It 
enables us to achieve almost perfect mode matching in 
the homodyne-detection, thereby making maximum use 
of the fragile non-Gaussian states. In this paper, we de- 
velop a theoretical time-domain description of the mode- 
filtering process, and then experimentally demonstrate 
the successful mode-reshaping as predicted from the the- 
ory. The simple technique demonstrated here can be 
straightforwardly applied to various photon-subtraction 
protocols, including state-generation jlo|, |lll| and entan- 
glement distillation [8]. 

Our mode-filtering strategy is especially effective for 
a wide range of CV protocols based on quantum tele- 



Noisy band 

Hi 



Bandwidth of 
operations 



ujQ—/\uj loq cjo+Aco' Freq. 
NT" 
Encoding 

(a) Sideband modes. 



.7 



^0 Freq. 

(b)A degenerate mode. 



c^o Freq. 

(c)A filtered degenerate mode. 



* Electronic address: takeda@ alice.t.u-tokyo.ac.jp| 
"I" Electronic address: ,akir af @ ap . t . u- 1 okyo .ac.jp] 



FIG. 1: (Color online) Frequency modes where quantum 
states are encoded. ujq\ the carrier frequency, 7: the band- 
width of the degenerate mode, k\ the cutoff of the optical 
high-pass- filter. 
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portation [121, [i3|, where the low frequencies are often 
contaminated by noisy measurements and feed forward 
operations |l3-[l6j. We have used the quantum telepor- 
tation apparatus in Ref. [7] as a testing device to show 
the advantages of mode-filtering. With this device, we 
teleported non-classical non-Gaussian states generated 
by mode-filtering, and then measured the degree of non- 
classicality preserved in the state. The Wigner func- 
tion of the teleported state has a minimum negativity 
of W(0,0) = -0.033 ± 0.005, compared to W(0,0) = 
—0.018 ± 0.007 without mode-filtering, where the nega- 
tivity here is an indication of the state's non-classicality 
{h = 1). These results show a clear filtering-related en- 
hancement of the quality, thereby demonstrating that our 
scheme is more effective for the robust implementation of 
various teleportation-based protocols with non-Gaussian 
states. 

This paper is organized as follows. In Sec. [Til we first 
review the basics of photon-subtraction and then intro- 
duce the idea and theory of mode-filtering. In Sec. ITTTl 
we give details of our experimental setup. In Sec. lIV| the 
experimental results of photon-subtraction with mode- 
filtering are presented and compared with the theory. In 
Sec.|Vl we show the results of quantum teleportation and 
demonstrate the advantages of mode-filtering. Finally 
Sec. IVII concludes this paper. 



II. SCHEME OF MODE-FILTERING 

Photon-subtraction is a method to conditionally gen- 
erate non-Gaussian states by the assistance of a single 
photon detection [2]. As illustrated in Fig. [21 a small 
fraction of a squeezed vacuum generated by a weakly- 
pumped OPO is reflected via a tapping beam splitter 
(TBS) towards an avalanche photo-diode (APD). This 
beam is used as a trigger beam, whereas the transmit- 
ted beam is called the signal beam. Conditioned on an 
APD click, one photon is subtracted from its correspond- 
ing quantum mode on the signal beam, and the resulting 
photon-subtracted state shows non-classical features in 
its non- Gaussian Wigner function with a negative dip. 

Since this scheme relies on the quantum correlation 
between the signal and trigger modes, as well as the pro- 
jective measurement of the trigger mode, the frequency 
mode of the non-Gaussian states is strongly correlated 
with the frequency mode of the photons detected by the 
APD. The degenerate OPO produces correlated photon 
pairs within a Lorentzian spectrum with bandwidth 7 in 
the weak pumping limit. Any photon within this spec- 
trum can be reflected by the TBS, trigger the APD click, 
and thereby herald a conditional photon in the signal 
mode. Thus, the correlation over the full spectrum is 
used to induce non- Gaussian states. As is proven in 
Ref. [17], the non-Gaussian states have the same fre- 
quency spectrum as that of the OPO, and their tem- 
poral mode function is given by the Fourier transform of 
the Lorentzian spectrum as fo{t) = -^e"^'*' (normalized 
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FIG. 2: (Color online) A simple schematic of photon- 
subtraction with quantum mode- filtering. APD: Avalanche 
photo-diode, CW: Continuous wave, HPF: High pass fil- 
ter, OPO: Optical parametric oscillator, TBS: Tapping beam 
splitter. 



as /|/o(t)|^(it = 1). In the time domain, this mode is 
a short wave packet of light, whose information can be 
extracted by multiplying /o(^) with the homodyne mea- 
surement signal. 

In this scheme, the frequency mode of the non- 
Gaussian state is totally determined by the Lorentzian 
spectrum of the OPO. From an engineering point of view, 
however, this frequency mode can be contaminated by 
low-frequency laser noise and also electric noise around 
the carrier frequency cjq (more details of the noise con- 
tamination are mentioned in Sec. IIIip . These types of 
noise are critical, and strongly deteriorate the fragile non- 
classicality of the non- Gaussian states. In order to avoid 
the noise contamination, we introduce a reflecting cav- 
ity before the APD. The cavity has bandwidth k which 
transmits photons in the low-frequencies and reflects the 
rest of the components towards the APD. The input- 
output relation of this cavity can be given in the fre- 
quency domain by [ISJ 

ar(cJo + cj) = — ain(^o + cl;) H — av(cJo + 

K — lUJ K — lUJ 



(1) 



where ain and are annihilation operators of the in- 
put and reflected modes, and is an auxiliary vacuum 
mode which does not affect the APD detection. The re- 
sponse function iuj/^n — iuj) in Eq. ([T]) shows that the 
reflecting cavity acts like an optical HPF^ which removes 
low frequency photons before photon-detection. As a 
result, quantum correlations in the low frequencies are 
unaffected by the projective measurement, and thus the 
non-Gaussian state is generated within a mode without 
low frequency components. 

Note that the optical HPF is introduced not in the 
signal mode to directly remove the low-frequency com- 
ponents of the state, but in the trigger mode to prevent 
the low-frequency quantum correlations from being used 
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for inducing non- Gaussian states. The generation rate 
of the non-Gaussian state is reduced as a result, but it 
does not directly affect the target quantum state itself: 
as we show below, the ideal heralded state is still the 
same non-Gaussian state. Therefore, this scheme makes 
it possible to reshape the frequency mode of the gener- 
ated state without any direct loss on or manipulation of 
the state within the mode. 

The mode-filtering scheme can be modeled in time- 
domain descriptions with an annihilation operator d{t) = 
(27r)~-^/^ / dujd{uJo + cJ)e~*^^ To begin with, a squeezed 
state generated from the OPO can be written in the time 
domain as 



exp(Pc - P'c) |0) 



n=0 



n! 



p^no), (2) 



where 



-i= J dtdt'C{t,t')a){t)a){t') 



(3) 



is the photon pair creation operator [19]. Cit^t') de- 
notes the two-time correlation function, which satisfies 
C{t,t') = C{t',t) = C{t - and therefore depends 
only on the time difference \t — because the output 
field of the OPO is stationary. Cit^t') can be shown to 
be 0(e), where e represents the pumping power of the 
OPO [l9|,|20|. 

We start with the simplest case of a weak pumping 
limit of e ^ 0. The higher order terms of n > 2 in Eq. ([2]) 
can be neglected to order 0(e). The squeezed state thus 
can be written approximately as |0) + (Pc'— ^c) 1^) ~ 
P(j) |0), and P(j |0) is the only term relevant to the photon 
detection. By the beam splitter transformation of the 
TBS a){t) VT^^dlit) + VRdl{t), |0) transforms 
into 

j dtdt'c{t,t') [(1 - R)dl{t)dl{t') 

+ 2^R{l-R)al{t)al{t') + Ral{t)a\{t')\ |0), |0), , (4) 

where the subscripts "s" and "t" denote the signal and 
trigger modes, respectively. Assuming P <C 1, we can ne- 
glect the case when more than two photons are reflected 
towards photon-detection. Thus, the only term in Eq. (|4]) 
relevant to the photon-detection has the form 



dtdt'C{t,t')dl{t)a\{t')\Q)^\Q), 



(5) 



Assuming an infinite bandwidth of the APD, we can 
model an APD detection at a time t = tc as a pro- 
jection onto (0|^at(tc)- This projective measurement 
of Eq. ^ gives / dtC{t - tc,0)at(t) |0)g. The operator 
can be regarded as a creation oper- 
ator of a single photon, whose temporal mode is defined 
by a mode function f{t) ex C{t — tcO)- In a sense, the 
beam splitter transformation and the projective measure- 
ment subtract one photon from P^t |0) by replacing the 



P(j operator with the corresponding single photon cre- 
ation operator. Thus, we can use fit) to observe the 
induced single photon state in the weak pumping regime 
of the conventional photon-subtraction. 

When the optical HPF of Eq. ^ is applied to the 
trigger mode before photon-detection, the effect can be 
described in the Schrodinger picture as the replacement 
oid\{(jjQ^uo) g{uj)dl{uJo-\-uj)^ where g{uj) = iuo / {n—iuj). 
In Eq. dnj, this effect appears as a convolution of C{t^t') 
with g{t) = (27r)-i/2 J duog{uo)e-''^^ as 



j dtdt' (^j drCit^t' 



Here the projection onto (0|^ ^t(tc) results in 



j dt(^j drCit- 



r)^(r))at(t)at(0|0)JO),. (6) 



(7) 



tc-r,0)^(-r) at(t)|0) 



In this case, the mode function can be defined as 
/HPF(t) (X j dTC{t-t,-T,())g{-T) cc fdrf{t-r)g{-r). 
Thus, the effect of the optical HPF can be included as 
a convolution of g{—t) with the original mode function 

fit)- 

A similar analysis can be applied to the higher or- 
der terms (n > 2) in Eq. (j2j). First, a simple induc- 
tive proof shows the following statement (A); "For any 
order n, (Pc - P^)"^ |0) can be decomposed into the 
terms written only with n or less P^ operators with var- 
ious correlation functions, including the vacuum term 
|0)." This statement (A) holds for n = explicitly. 
Suppose that the statement (A) holds for n = m, and 
therefore (Pc - P^t)^ |0) can be decomposed into the 
sum of the terms in the form |0) or P^t^P^^ • • • P^^ |0) 
{k < m), where Ci{t^t') {i = l,2,...,/c) is the cor- 
relation function of each operator. In this case, all 
of the terms resulting from the (m + l)th order term 
{Pc-Pc)'^^^ |0) can be obtained by applying (Pc-P^) 
to each mth order term as {Pc — P^t) |0) = —Pq |0) oi" 

(Pc - Pc)PcrPc2 '"^L 1^)- ^^^^^^ 

composed with the commutation relation [d{t)^ d^f)] = 

5(t-t') and [d{t),d{t')] = [at(t),at(tO] = as 



iPc-Pl)Pl,Pl,---P'cjO) 

k 



^ i 
k 

4 V At At •••At P^ •••P^ P^ •••P^ 

i<j 

-P'cPk---Pk)\^)^ (8) 



where q = j dtdt' C*{t,t')Ci{t,t') and Cij{t,f) = 
f dTdT'C*{T,T')Ci(t,T)Cj(T' ,t'). Here, the second term 
in the right-hand side of Eq. (|8]) vanishes when k = 1. 
These decompositions show that all of the (m + l)th 
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order terms can also be decomposed into the form 



pt pt 



• P^, |0) (/c < m + 1), and thus the statement 



(A) holds for n = m + 1 as well. The above discussions 
prove that the statement (A) holds for any order n. 

Then, in a similar way, the transformation of the TBS 
and the projection onto (0|^at(tc) subtract one pho- 
ton from one of the correlated photon pairs in each 
term P^t^ • • • P^t^ |0). As a result, various photon- 
subtracted terms are produced by the replacement of 
one of these P^ operators with its corresponding sin- 
gle photon creation operator. The replacement of P^. 
(i = 1,2, gives 



j dtCi{t-tc,0)al{t) 



Here, the projective measurement induces a single pho- 
ton defined by the mode function fi{t) ex Ci{t — tc^O). By 
selecting the fi{t) mode, we can observe the pure single 
photon induced by the APD detection, as well as some 
fraction of the photon pairs generated in the OPO. When 
applying the optical HPF, we need to replace Ci{t — tc^ 0) 
in Eq. 1^ with / drCi{t - - r,0)^(-r). The effect of 
mode-filtering here can also be included as a convolution 
of fi{t) with g{—t), as shown previously. Note that dif- 
ferent terms resulting from the projection can have sin- 
gle photons with different mode functions (when i 7^ j, 
fi{t) fjit) in general), thereby making it impossible to 
extract all the induced single photons purely at the same 
time. However, the mode-filtering effect on all the mode 
functions /2(t ),..., //c(t)} can be written in the 

same way as the convolution with g{—t). Therefore, when 
we introduce the optical HPF to the photon-subtraction 
scheme, we only need to replace the originally-used mode 
function f{t) with J drf{t — r)g{—r), regardless of the 
pumping power e. 

Our experiment is performed in the weak pumping 
regime, thus we can assume that the original mode func- 
tion is fo{t) = y^e~^l^l. g{t) can be derived from the 
Fourier transform of g{uj) = iuo / — iuo) as 



g{t) V27r X 




m) {t>0) 
{t<0) 



(10) 



Thus, by convoluting fo{t) with g{—t)^ and then set- 
ting the normalization constant properly, we obtain the 
optically-high-pass-filtered mode function as 



/HPF(t) = X 




2k ^Kt 
7 — K 



(t>0) 

(t < 0) ■ 



(11) 



The temporal shapes of Eq. (pTj) are shown in Fig. [3l 
Contrary to the original non-negative mode function, this 
function has a negative dip due to its lack of low fre- 
quency components. 
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FIG. 3: The dependence of the optically high pass filtered 
mode function fuPF{t) on cutoff k. We set 7 = 27r x 6.2MHz. 



III. EXPERIMENTAL SETUP 

Our experimental setup consists of three parts; the first 
part is for photon-subtraction, the second part for quan- 
tum teleportation, and the third part for evaluating the 
quality of the teleported states (Fig. S]). Apart from the 
optical HPF part, the setup is mostly identical to that in 
Ref. [7], where more details are provided. 

The output of a single-mode continuous-wave 
Tiisapphire laser (SolsTiS-SRX, M Squared Lasers) 
around 860 nm is used as the light source, a part of 
which is frequency-doubled for pumping three optical 
parametric oscillators (OPO). First, in the photon- 
subtraction part, a weakly squeezed vacuum is generated 
from OPOl (half width at half maximum: 6.2 MHz) 
with a 15 mW pump. The output of OPOl is directed 
into a TBS which reflects 3% of the beam toward an 
APD as the trigger beam, while the other beam is 
sent to the teleportation setup as the signal beam. 
In the trigger channel, photons in all of the modes 
except for the degenerate squeezed mode are blocked 
with two Fabry-Perot filtering cavities with much wider 
HWHM (55 and 18 MHz) than that of OPOl. Then 
the trigger beam is injected into an optical HPF, which 
is a Fabry-Perot cavity with a round-trip length of 
280 mm (HWHM: 500 kHz, Finesse: 1080). Here the 
theoretical cutoff k is set to 27r x 0.5 MHz. By increasing 

we can reduce the low-frequency noise contamination. 
However, this decreases the event-to-dark-count ratio on 
the APD, which degrades the purity of the generated 
states. Therefore the choice of is a compromise 
between these two effects. The reflection of the optical 
HPF is mode-reshaped and sent to the APD, while the 
transmission is monitored on a photodiode for locking 
the cavity. Ideally, the cavity perfectly blocks the ujq 
component when locked on resonance, but in our setup 
14% of the reflected component remains due to a slight 
mismatch between the transmittances of two cavity 
mirrors (impedance mismatch). This effect can be taken 
into account by a more detailed model of the mode func- 
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FIG. 4: (Color online) Experimental setup. APD: Avalanche 
photo-diode, EOM: Electro-optic modulator, HBS: Half beam 
splitter, LO: Local oscillator, OPO: Optical parametric oscil- 
lator, PBS: Polarizing beam splitter, QWP: Quarter wave 
plate. 



nal for the evaluation. These types of noise often need 
to be removed with electrical HPFs. In our setup, a 
HPF with a cutoff frequency of 101 kHz is required in 
the final homodyne signal. Non-Gaussian states gener- 
ated without mode-filtering have low-frequency compo- 
nents, which are inevitably removed by this HPF. The 
resulting frequency mode mismatch in the homodyne- 
detection directly deteriorates the non-Gaussian states, 
thereby limiting the quality of the teleportation. In con- 
trast, non-Gaussian states generated with mode-filtering 
do not have frequency components lower than 101 kHz, 
thus we can avoid the low-frequency noise contamination, 
as well as achieve almost perfect mode matching in the 
final homodyne-detection. 



IV. NON-GAUSSIAN STATE GENERATION 
WITH MODE-FILTERING 



tion, which is described in the Appendix. Conditioned 
on an APD click, a photon-subtracted squeezed vacuum 
is generated in the signal beam and used as an input 
state for the subsequent quantum teleportation circuit. 

For realizing quantum teleportation of this state, first 
Alice and Bob need to share a pair of broadband EPR- 
correlated beams as a resource. This pair can be obtained 
by mixing two strongly squeezed vacua from 0P02 and 
0P03 (HWHM: 12 MHz) on a half beam splitter. Alice 
combines the input state with one of the EPR pairs on 
a half beam splitter, then measures the x and p quadra- 
tures of the resulting two outputs. The measurement 
results are sent to Bob as gain-tuned electric signals, 
which Bob uses to modulate displacement beams (shown 
in Fig. HI and to correct the other beam of the EPR pair. 
As a result, the input state is reconstructed. In order 
to faithfully teleport the quantum state within the origi- 
nal wave packet quantum mode, all these operations are 
performed over a frequency spectrum broader than this 
wave packet spectrum. 

In order to evaluate the quality of the teleported quan- 
tum states, the output beam is measured by a balanced 
homodyne detector. The local oscillator phase in this de- 
tector is scanned for the full angular resolved tomography 
of the state, and 100,000 quadrature data are recorded 
over [0,27r]. We use the inverse Radon transform based 
on polynomial series expansion [2l] to reconstruct the 
Wigner functions W{x^p). 

In this teleportation circuit, and in other teleport ation- 
based circuits as well, low-frequencies can be contami- 
nated by mainly two sources of noise. The first source 
is the laser noise around the carrier frequency. The non- 
Gaussian states can be deteriorated by this noise derived 
from the local oscillator beams in imperfectly-balanced 
homodyne detections for Alice and the final evaluation, 
as well as the displacement beams directly mixed by 
Bob. The second source is the low- frequency electric 
noise, which contaminates Alice's homodyne signal for 
the displacement operation and the final homodyne sig- 



First the states generated by photon-subtraction are 
analyzed to demonstrate the mode-filtering effect by the 
optical HPF. We record 200,000 frames of quadrature 
data both in the conventional and mode-filtering photon- 
subtraction. When using the conventional method, we 
unlock the optical HPF and monitor the signal of the 
transmitted beam to keep it totally off resonant during 
the measurement time. 

Experimental Wigner functions W{x,p) in these two 
methods are shown in Fig. |5(a)| and |5(b)[ Here we 
applied the mode function including the effect of the 
impedance mismatch of the optical HPF (Eq. ([H]) in 
the Appendix). Each Wigner function is reconstructed 
by optimizing the parameters 7 and of the mode func- 
tion to extract the minimum 1^(0, 0) {hi is set to for 
the conventional method). The optimal parameters of 
7 = 27r X 6.5 MHz and = 27r x 0.48 MHz match well 
the experimentally estimated values of 27r x 6.2 MHz 
(OPOl bandwidth) and 27r x 0.50 MHz (optical HPF 
bandwidth). Both Wigner functions have minimum neg- 
ativities of W(0, 0) = -0.171 ± 0.003 (conventional) and 
W(0,0) = -0.179 ± 0.003 (mode-filtering) without any 
corrections, thereby excluding any description as a clas- 
sical probability density. The experimental Wigner func- 
tions sometimes have minimum values not exactly at the 
origin, but at a slightly shifted point around the origin 
possibly due to the interference of the trigger photons 
with the other unwanted fake trigger photons. Here, 
we adopt the minimum negativity around the origin as 
W(0,0). This non-classical negative region is fragile and 
easily disappears with every percent of experimental loss 
and other imperfections. Thus, the quality of these non- 
Gaussian states can be evaluated by their negativity. 

In Fig. |5(c)[ we plot the dependence of the minimum 
value VK(0, 0) on the cutoff of applied mode functions 
(7 is fixed at 7 = 27r x 6.5 MHz). Without the opti- 
cal HPF, the theoretical mode function fo{t) = y^e"^'^' 
corresponds to the case of /t: = in Eq. (p!4|) . There- 
fore W(0, 0) increases with increasing as a result of the 
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(a)A Wigner function with the 
conventional method. 7/27r = 6.5 MHz, 
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(c)The dependence of W(0, 0) on the cutoff k of 
the applied mode function (7/27r = 6.5 MHz). 
Theoretical curves based on Eq. (|12|) are plotted 
as well. 



FIG. 5: Experimental results of photon-subtraction with or without mode- filtering. 



mode mismatch between the photon-subtracted temporal 
mode and the applied mode function. With the optical 
HPF, the plot shows the different hz dependence which 
has an optimal K:. 

These experimental data are in good agreement with 
the theoretical curves calculated from a simple and realis- 
tic model [22I , thereby showing the experimental success 
of the mode-filtering (Fig. |5(c)' ). In the theoretical model, 
the generated state is modeled by aS{r) |0), where a is 
the annihilation operator describing photon subtraction 
and S{r) is the squeezing operator with a squeezing pa- 
rameter r. This state is degraded by experimental losses 
and imperfections, and as a result the value of the Wigner 
function at the origin can be expressed as |22j 



1 - 27^ + 2Cr] [1 + 2(1 - T]) sinh^(7 



TT [1 + 47^(1 - T]) sinh^(r) 



1 3/2 



(12) 



where rj is the overall efficiency and ( is given as a dark 
count rate divided by a total count rate. The ideal limit 
of 7^ = 1 and C = gives W(0, 0) = — I/tt regardless of r. 

From our analysis the squeezing parameters were es- 
timated as r = 0.38 ± 0.02 (conventional) and r = 
0.36 ± 0.02 (mode- filtering). The total count rates were 
8200 ± 300 and 7200 ± 300 per second respectively, in 
which the dark count rate of 150 ± 30 per second was 
included. Therefore ( = 0.02 for each case. The over- 
all efficiency can be calculated as 7/ = r]o{r]^)'^^ where 7^0 
is the efficiency before applying a mode function, and 
Vk. = J fopt{t)fK{t)dt is a mode matching parameter be- 
tween the optimal mode function /opt(^) and the applied 
mode function /,^(t). Using the parameters above, we 
calculated VK(0, 0) for each by fixing the unknown pa- 
rameter 7^0 for the most plausible value of 0.83. 

7^0 was also calculated directly from experimentally es- 
timated values as 7/0 = ^op0^pr(^vis)^^hom, where the 
OPO escape efficiency 7^opo = 0.98, the propagation ef- 



ficiency r]pr = 0.95, the visibility 7^vis = 0.98, the homo- 
dyne efficiency 7/hom = 0.95 which includes the efficiency 
of the photodiode and the electrical signal-to-noise ratio. 
All these values gave r]o =0.85, which is reasonably close 
to the directly estimated value of 0.83. 

It appears there is some discrepancy between experi- 
ment and theory at z^: = in Fig. |5(c)[ which we partly 
ascribe to some low frequency noise. It implies that the 
negativity in the conventional method is partly limited 
by this noise, whereas in the mode-filtering method the 
mode does not have low frequency components and thus 
can avoid this limitation. 

These data and analysis show the successful mode- 
filtering as predicted from the theory. Therefore, we suc- 
ceeded in tailoring the quantum mode of non-Gaussian 
states, while preserving or even enhancing the non- 
classicality of the state within the mode. 



QUANTUM TELEPORTATION WITH 
MODE-FILTERING 



The next step is to demonstrate the advantage of 
our mode-filtering method by measuring the degree of 
non-classicality preserved in the state after teleporta- 
tion. For this purpose a tomography data set of the 
teleported states with and without the mode-filtering 
(100,000 frames each) is measured 4 times in a row. We 
reconstructed Wigner functions from each data set with 
optimized parameters of mode functions. The minimum 
negativities W(0, 0) of all data sets are shown in Fig. |6(a)[ 
and a selecti on of the teleported Wigner functions are 
shown in Fig. |6(b)| and |6(c) 



All of the Wigner functions exhibit negativities around 
the origin, showing the successful teleportation of the 
fragile non-classical feature. Negative values are es- 
timated from four sets of reconstructed Winger func- 
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(a)VK(O, 0) of four data sets. Averaged values are 
plotted with dashed (mode filtering) and dot-dashed 
(conventional) lines. 
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(b)A Wigner function with the 
conventional method. 
W(0, 0) = -0.017 ±0.005. 
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(c)A Wigner function with the 
mode-filtering method. 
1^(0,0) = -0.033 ±0.005. 



FIG. 6: Experimental results of quantum teleportation with or without mode- filtering, (b) and (c) correspond to "Data set 
number 2" of (a). 



tions as W(0, 0) = -0.018 ± 0.007 (conventional) and 
W(0,0) = -0.03 3 ± 005 (mode-filtering). The Wigner 
functions in Fig. |6(b)| and 6(c)| exhibit a clear difference 
in the depth at the center dip. The mode-filtered states 
exphcitly preserve stronger negativities, demonstrating a 
clear filtering-related enhancement of the states' robust- 
ness in the teleportation process. Due to the impedance 
mismatch of the optical HPF, our mode-filtered states 
still have some low- frequency components, and are there- 
fore partly affected by the noise contamination. If per- 
fect impedance matching is obtained, the enhancement 
can be even greater. 

In summary, we have experimentally demonstrated 
that our efficient encoding strategy with mode-filtering 
is highly effective in the current teleportation appara- 
tus, where low frequencies are contaminated by the laser 
noise and electric noise in three homodyne measurements 
and two displacement operations. This filtering-related 
enhancement can also be seen in a wide range of CV pro- 
tocols, especially teleportation-based protocols with the 
same kind of inevitable noise contamination in practice. 



VI. CONCLUSION 

In conclusion, we have developed and demonstrated 
a mode-filtering technique of no n- Gaussian states gener- 
ated by photon-subtraction. By introducing a reflecting 
cavity in the trigger channel, we removed the problem- 
atic noisy low- frequency components of the non- Gaussian 
states. The mode-reshaped states are analyzed with the 
mode function derived from a theoretical model, and the 
agreement of experiment and theory shows the success- 
ful mode-filtering without any degradation of the non- 
classical states within the mode. After the teleportation, 
the mode-filtered states show stronger non-classicality 
than the states without mode-filtering, thereby show- 



ing that our technique is effective to achieve robust im- 
plementation of CV protocols, especially teleportation- 
based protocols utilizing non-Gaussian states. This tech- 
nique can be straightforwardly applied to various types 
of photon-subtraction protocols, and thus it should be 
an important ingredient in a variety of applications of 
non- Gaussian states towards universal CV QIP. 



APPENDIX: THE EFFECT OF IMPEDANCE 
MISMATCH 

In general, the reflected mode of the optical HPF can 
be written in the frequency domain as [18|] 

. / . {i^i - ^2)/2 + icj ^ , , 



{k\ + ^2)/2 - iijJ 



ttv(^0 



(13) 



where k\ and K2 denote the decay rates of two cavity 
mirrors. If hci = f<i2 = the impedance match is perfect 
and Eq. (p!3|) gives Eq. ([T]). Otherwise we define n = 
{hzi + A>:2)/2 and hz^ = {hzi — /^2)/2, then the response 
function of the cavity is given by g{uj) = {n' ^iuo) / {n—iuj) . 
Since ^(0) 7^ 0, some part of the uoq component is reflected 
toward the APD. By convoluting g{—t) with the original 
mode function fo{t) = y^e~^'^', we obtain 



fiiPF{t)=Nx 



(t>0) 



where A/" is a normalization constant given by 



N ■ 



(14) 



(15) 
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In this experiment, k.' is set to k.' = — 0.37a^ based on and FIRST commissioned by the MEXT of Japan, and 
experimental values. ASCR-JSPS. The authors acknowledge Seiji Armstrong 



for his assistance with the editing of this paper. 
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